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Abstract 

Within the inflationary scenario, Planck scale physics should have affected the comoving modes' initial 
conditions and early evolution, thereby potentially affecting the inflationary predictions for the cosmic 
microwave background (CMB). This issue has been studied extensively on the basis of various models for 
how quantum field theory (QFT) is modified and finally breaks down towards the Planck scale. In one 
model, in particular, an ultraviolet cutoff was implemented into QFT through generalized uncertainty 
relations which have been motivated from general quantum gravity arguments and from string theory. 
Here, we improve upon prior rmmerical and semi-analytical results by presenting the exact mode solutions 
for both de Sitter and power-law inflation in this model. This provides an explicit map from the modes' 
initial conditions, which are presumably set by quantum gravity, to the modes' amplitudes at horizon 
crossing and thus to the inflationary predictions for the CMB. The solutions' particular behaviour close 
to the cutoff scale suggests unexpected possibilities for how the degrees of freedom of QFT emerge from 
the Planck scale. 
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1 Introduction 

Within the inflationary scenario, spacetime typically 
inflated to the extent that comoving modes which are 
of cosmological size today originated well beyond the 
Planck scale (and any other natural ultraviolet (UV) 
cutofl' scale such as the string scale). This observa- 
tion has led to a growing body of work exploring the 
tantalizing prospect that the predictions of inflation 
could be observably affected by quantum gravity. For 
early papers on this subject see, e.g., 

Technically, all possible Planck scale effects on 
the comoving modes, fc, at late times take the form 
of a nontrivial fc-dependent selection of solution to 
the ordinary mode equation, see, e.g., This is 
because, independent of the details of Planck scale 
physics, each mode's evolution equation reduces to 
the ordinary wave equation as soon as its proper wave 
length is significantly larger than the Planck scale. It 
has been possible, therefore, to determine on these 
general grounds that quantum gravity effects should 
manifest themselves qualitatively in the form of su- 
perimposed oscillations in the CMB spectra as well 
as in a breaking of the consistency relation between 
the scalar/tensor ratio and the tensor spectral index 
of single field inflation, see, e.g., |71 151 1^ ITUl [TT| . 

Any quantitative predictions, however, are par- 
ticular to each particular theory of quantum gravity. 
Concretely, for wavelengths closer and closer to the 
cutoff scale the framework of QFT should hold - with 
characteristic increasing corrections - until it flnally 
breaks down at the cutoff scale. A theory of quantum 
gravity, be it, e.g., loop quantum gravity or M theory, 
should eventually allow one to calculate those charac- 
teristic corrections to the framework of QFT in this 
"sub-Planckian" regime, i.e., in the regime of wave- 
lengths larger than but close to the cutoff length. In 
addition, the theory of quantum gravity should allow 
one to determine the initial condition for the evo- 
lution of a comoving mode within the framework of 
QFT, namely at the time when the mode's proper 
wavelength first exceeds the cutoff length. 

In this context, first-principle calculations are very 
difficult, of course. Therefore, a number of authors 
have modelled the corrections to the framework of 
QFT from quantum gravity simply as UV modifica- 
tions to the dispersion relations. The implied effects 
on the inflationary predictions for the CMB have been 
investigated, for example, in ^EIEIIEI- Particu- 
lar modiflcations to the dispersion relations for short 
wavelengths have been motivated either ad hoc or on 
the basis of analogies to the propagation of waves in 
condensed matter systems. It was found that the ef- 
fects of Planck scale physics could be considerable. 



even though there appear to be strong constraints 
due to a possibly strong backreaction problem, see, 
e.g., d. 

A less ad hoc approach 0] has been to model 
the behaviour of quantum field theory in the sub- 
Planckian regime through the implementation of cor- 
rections to the first quantization uncertainty rela- 
tions, see 25 . These corrections implement an ultra- 
violet cutoff in the form of a finite minimum uncer- 
tainty in spatial distances, which has been motivated 
by studies in quantum gravity and string theory, see, 

e.g., [isiiiaiiiiiiiiiini. 

This approach to modelling Planck scale physics 
in inflation was further investigated in particular in 
[3 [3 El d] • So far it has been possible to solve only 
approximations of the mode equation that arises in 
this approach. In fact, to solve the exact mode equa- 
tion is known to be highly nontrivial even numeri- 
cally because of a particular singular behaviour of the 
mode equation at the mode's starting time. Here, we 
are following up on this series of papers. 

In particular, we calculate here the exact solu- 
tions to the exact mode equation in both the de Sit- 
ter and power-law backgrounds. This allows us then 
to follow the mode solutions back towards the time 
when the proper wavelength approaches the cutoff 
length. That in turn allows us to calculate the ex- 
act behaviour of physical quantities such as the field 
fluctuation amplitudes and the Hamiltonian towards 
the Planckian regime. We obtain, therefore, a fully 
explicit model in which to explore possible mecha- 
nisms that could impose initial conditions on comov- 
ing modes as they enter into the regime of validity of 
quantum field theory. 

2 Short distance physics and in- 
flation 

2.1 UV cutoff through minimum length 
uncertainty 

A rather general assumption about quantum field 
theory in the sub-Planckian regime is that it should 
still possess for each coordinate a linear operator x* 
(inherited from first quantization) whose formal ex- 
pectation values (e.g., in the space of fields that are 
being summed over in the path integral) are real. The 
X* may or may not commute. As shown in 14^, this 
implies that the short distance structure of any such 
coordinate, considered separately, can only be contin- 
uous, discrete, or "unsharp" in one of two particular 
ways. Studies in quantum gravity and string the- 
ory point towards one of these two unsharp cases, 
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namely the case of coordinates x whose formal un- 
certainty Ax possesses a finite lower bound Axmin at 
the Planck or string scale. This short distance struc- 
ture arises from quantum gravity correction terms to 
the uncertainty relation (see |2(J)'): 

AxAp>^{l + f3{Apf + ...) (1) 

Here, /3 is a positive constant which, as is easily veri- 
fied, parametrizes the cutoff length through Axmin = 
V^. The units are such that h = 1. As was first 
pointed out in j2f| , uncertainty relations of this form 
arise from corrections to the canonical commutation 
relations which, for example in the case of one space 
dimension, take this form: 

[x,p]-i(l + V + (2) 

It was shown in j22| that the multi-dimensional gen- 
eralization is unique to first order in P if rotational 
and translational symmetry is to be preserved. 

Let us remark that, as was shown in |23| . this type 
of natural UV cutoff implies that physical fields pos- 
sess a finite bandwidth in the information theoretic 
sense: In the presence of this short distance structure 
one can apply Shannon's sampling theorem to con- 
clude that if the values of a physical field are known 
on any set of points with average density above twice 
the Planck density then the values of the field every- 
where are already determined and can be explicitly 
calculated from these samples"'^. Notice that no sam- 
pling lattice is preferred, as all sampling lattices of 
sufficiently high average density can be used to recon- 
struct the fields. Therefore, in a theory with this type 
of natural cutoff spacetime can be viewed as continu- 
ous, in which case the conservation of the spatial sym- 
metries is displayed, while, fully equivalently, in the 
same theory spacetime can also be viewed as discrete, 
in which case the theory's UV finiteness is displayed. 
So far, this type of natural UV cutoff has been con- 
sidered mostly within models that break Lorentz in- 
variance, as we will here. The generalization to a 
generally covariant setting has been started in |24j. 

2.2 Mode generation in the presence 
of the UV cutoff 

In 01 , the short distance cutoff formulated in Q was 
applied to the theory of a minimally coupled massless 

'^The situation is mathematically identical, for example, to 
that of HiFi music signals of bandwidth 20KHz, whose ampli- 
tude samples are recorded on CDs at a rate of 4 ■ 10'' samples 
per second. Using Shannon's theorem, CD players are able to 
precisely reconstruct the continuous music signal at all time 
from these discrete samples. 



real scalar field (/>(x, t) in an expanding Friedmann- 
Robertson- Walker (FRW) background. This setup 
can be used to describe the quantum dynamics of the 
tensor as well as the scalar fluctuations in inflation. 
Concerning the mode evolution, the only difference 
is that, in the latter case, the scale factor a in the 
mode equation is replaced by z = (/)ga^/a'. Here, 4>o 
is the bulk zero mode of the inflaton field with the 
prime denoting differentiation with respect to ij. (It 
should be kept in mind, however, that the Hamilto- 
nian of the gravitational waves that are due to tensor 
modes contributes to the energy momentum tensor 
only from second order.) 

The approach of @J was continued in [S] as well 
as in a series of papers by Easther et al. [71 |H1 E| • In 
each case, for simplicity, the notation for the tensor 
modes was used (i.e. with the mode equation con- 
taining a instead of z), as we will also do here. 

We assume the case of spatial flatness of the back- 
ground spacetime and we use the conformal time co- 
ordinate -q. Note that the UV cutoff Axmin = \fP is 
introduced in proper distances as opposed to comov- 
ing distances. The action for the field (\>{x, rj) then 
reads, see 0]: 




Note that our function 4>f.{ri) is related to in [7] 
by uj^ = Xja} ■ 0j,(?/) and that our function 1/(77, ^) is 
the function v(r]. p) from 7 multiplied by a'*. Using 
the operators 

x' : ^{x,r\) Xi(t){x,r]) 

p' : 0(a;,ry) -id/dxi(j}{x,r]) 

the action can be expressed representation indepen- 
dently in terms of these operators and the usual 
inner product of fields. While the fields' expression 
for the action ||3Jl is held fixed the underlying com- 
mutation relations are now corrected to introduce the 
desired cutoff Q . The modification of the short dis- 
tance behaviour then manifests itself whenever one 
writes the action in any representation, such as in 
comoving momentum space. For example, at dis- 
tances close to the cutoff length, comoving modes 
no longer decouple. However, it is then still possi- 
ble to find variables (77, k) in which the k modes of 
the fields' Fourier transform do decouple^. The de- 

•^Note that, as discussed in and these variables k 
approximate the comoving momentum variables k only at long 
wavelengths. 
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coupling modes are important as they describe the 
truly independent degrees of freedom. They will be 
denoted as 4>j.{'i]). The realness of the field (j){'K,ri) 
translates into 0* — (f)_^. The action |j3Jl can then be 
written in the form 



with the Lagrangian 



where the coefficient functions are given by 

a^plog(— /3fc^/a^) 



v{ri,k) 



(3 l + plog(-/3fc2/a2) 



exp 



blog(-/3P/a') 



1 + plog(-/3A:2/a2 



(4) 



(5) 



(6) 



(7) 



The function plog used here, the "product log", is 
defined as the inverse of the function x xe^ (also 
known as the Lambert W-function). It will be impor- 
tant later that, at x = — 1/e, this function possesses 
an essential singularity. 

Equally important will be the way in which the 
short distance cutoff^ affects the integration region 
of the action functional Q): each mode fc enters the 
action once a(r\) has grown enough so that the con- 
dition fc2 < c? jefi holds, i.e., when: 

a(?7c) = k\fe^ k, kAxjnm- 

Correspondingly, the mode equation that follows from 
the action 0, see 0], 



-I- 
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^(8) 

is an equation of motion which possesses a starting 
time. As had to be expected, the implementation of a 
UV cutoff in an expanding spacetime has led us to the 
problem of setting initial conditions on these newly 
emerging comoving modes. Before we solve the mode 
equation to obtain the solutions' early time behavior, 
let us briefly consider the mode solutions' late time 
behaviour. 



^For /3 — > 0, 13 reduces to the standard action. Note that 
Fourier transforming and scahng only commute up to a scahng 
factor of a^. 



2.3 Late time behaviour 

At late times, when a mode's proper wavelength far 
exceeds the cutoff length, the mode becomes insen- 
sitive to the nontrivial short distance structure of 
spacetime. The modified mode equation © then re- 
duces to the mode equation without cutoff, which 
(taking into account the non-commutativity of scal- 
ing and Fourier transforming) reads with our conven- 
tions'*: 

Let us derive the late time behaviour of namely 
^ 0, in de Sitter space a{r]) = ^l/Hrj. Since in 
this case we have — = —- and — — while k^ is 
constant, equation Q on large scales becomes: 



4I, 



= 



The solution to this equation is given by 



B2 

„3 " 



(10) 



So for late times, 77 — > 0, the mode function 4>j^{ri) 
diverges cx ^ in the de Sitter case. We recall that 

= o.^{v)^k{''l)- Thus, as expected, the physical 
field's mode function $^(77) approaches a constant for 
late times. Let us remark that the analogous calcu- 
lation for classical fluctuations indicates that classi- 
cal scalar or tensor fluctuations that existed before 
inflation would have their wavelength stretched but 
their amplitude maintained, i.e., they would not be 
"ironed out" by inflation. Inflation is usually thought 
to have lasted long enough so that, classically, the 
inflation-induced flattening of a fixed proper spatial 
volume would hinge on the assumption that there 
were no pre-existing spatial "ripples" down to wave- 
lengths shorter than the Planck length. Of course, 
the Planck scale is expected to rule out the notion 
of pre-existing classical ripples that are shorter than 
the Planck length. 



3 The Hamiltonian 

Using the canonically conjugate field to see 

'^kiv) = '^4'' ^iv) ^ '^'^^4'-k{v)i the Hamiltonian cor- 
responding to the action takes the form 



1 



1 



(11) 

^Recall that at large scales we also have k ^ k. 
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To quantize one promotes (fif, and tt^ to operators 
satisfying the usual commutation relation [(f)f,,TTf] = 



in momentum space. Note that the field 



commutators are of course modified when expressed 
in position space. 

In the Hamiltonian, the integration region increa- 
ses over time, expressing that the field operators of 
comoving modes enter the Hamiltonian only when the 
mode's wavelength exceeds the cutoff length. This 
is significant because it means that the Hamiltonian 
contains operators which act nontrivially on the Hil- 
bert space dimensions of a particular comoving mode 
k only from that time rjdk) when this mode's proper 
wavelength starts to exceed the cutoff length. In the 
literature, this time (which coincides with the essen- 
tial singularity of the plog function) has been called 
the mode's "creation time" and we will stay with 
this terminology. Let us keep in mind, however, that 
the dimensions on which the harmonic oscillator of 
a given mode is represented are of course always in 
the Hilbert space. In the Heisenberg picture, all that 
happens at the creation time of a mode is that the 
mode's field and conjugate field operators enter the 
Hamiltonian for the first time, while conversely, in 
a shrinking universe, a mode's field operators would 
drop from the Hamiltonian when the mode's wave- 
length drops below the cutoff length. 

We express the quantum field in terms of creation 
and annihilation operators at and and in terms of 



the mode functions 0^(?7): 



(12) 



Note that the notation 4>i^{ri) now no longer stands 
for the classical field (which obeys 0| = (p^k) but in- 
stead for a mode function. For the mode functions we 
can choose, as usual that (p^ — The Wronskian 
condition 



,^{v,k) [MvWiiv) - <l>l{vWv)] = (13) 

insures that the so-defined quantum fields obey the 
canonical commutation relations. We can then re- 
arrange the quantized Hamiltonian in the following 



form that we will later refer to 
H = 



2 



X («fca-fc) 



5\Q). (14) 



To see that this Hamiltonian gives the correct 
equation of motion JHJ, it is straightforward to ap- 
ply the Heisenberg equations and use the definition 
of TT^. Note again the bound on the integral, which 
reflects the fact that each mode k only contributes 
to the Hamiltonian from its creation time rjdk) on- 
wards; before that time, the operators and at are 
not contained in H . 

Let us note that, curiously, in the case of a shrink- 
ing universe the modes' creation and annihilation op- 
erators , at successively drop out of the Hamilto- 
nian. Thus, a growing number of states can no longer 
be "reached" by the Hamiltonian and the field oper- 
ators while the universe is contracting. All informa- 
tion encoded in these states would no longer inter- 
act although it would become relevant again during 
a subsequent expanding phase. 

We will be interested, in particular, in the ground 
state energy term of each comoving mode, given by 
the last line of 1(1^ : 



:.(») 



(15) 



\4>kiv)\ m(»7)-9 



a{T]) 



Since the Hamiltonian is closely related to the 00 
component of the energy momentum tensor T'^^, (|15() 
essentially describes the energy contribution of each 
mode which arises when the mode outgrows the cutoff 
scale. Two time dependencies determine the evolu- 
tion of the modes' ground state energy (|15|) : first, 
the behaviour of the coefficients v{ri,k) and ii{r],k), 
and second the mode function 4'j.{ri). While the for- 
mer are already known explicitly (see JHl and (|7|)). as 
we will discuss, the choice of solution to the equation 
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of motion |(SJ) for the mode function will nontrivially 
affect the ground state energy. 

4 Solving the mode equation 

We will now derive the exact solutions to JSJ which 
obey the Wronskian condition (|13|l . The difhculty 
here is due to the plog-functions in ii{ri, k) and i'{'q, k): 
the equation of motion for each fc-mode has an irreg- 
ular singular point at the mode's creation time r]c{k). 
In previous work, [7] , the mode equation was replaced 
by an approximate mode equation which had a reg- 
ular singular point at the mode's creation time. The 
exact solutions to this approximate mode equation 
were obtained both as power series and in analytical 
form. 

Our aim now is to obtain the exact solutions to 
the exact mode equation. To this end, we employ 
a suitable variable transformation that turns the ir- 
regular singular point into a regular singular point, 
without the need for approximations. We then apply 
the Frobenius method to derive the exact solutions, 
namely in terms of convergent power series and log- 
arithms. 

4.1 Singularity transformation 

Technically, ric{k) is called an irregular singular point 
of ijHI) because the term 



and the term 



{v - Vcik)) ■ — 



do not have convergent Taylor expansions around this 
point. Let us now define a new variable r. 



(16) 



which plays the role of time, but is a function of both 
f] and k. The scale factor 0(77) and the coefficient 



functions 1^(77, k) and ^{rj, k) then read^ 
air) 
fi{T,k) := 
i>(T,fc) := 



T ' 

e-(l + r)2 



/32fc4(l-Hr) 



(17) 
(18) 
(19) 



It is straightforward to translate the derivatives d/ drj 
into derivatives in terms of derivatives d/dT and we 
then obtain that the equation of motion ((HJ takes the 
form: 



= 





v"\ 










^.1^ - 3 




)- 








au 


a 





(20) 



Note that we now let the prime denote differentiation 
with respect to r rather than 77. 

In order to make those expressions explicit in an 
example, let us consider the case of de Sitter space 
a (77) = —l/Hrj. In this case, by inverting (|16|l . 7/ is 
expressed in terms of r and k through: 



77(r, k) 



1 V=TC^ 
Hk VP 



(21) 



Using lfT7|) . (fTH|l and (fT^ . we then obtain for the 
mode equation in de Sitter space: 



1 + 15g2/3 + 9H^(3t 



2r(H-r) 



(22) 

The creation "time" (when a mode outgrows the min- 
imum length) now corresponds for all fc to Tc = —1. 
(This is particular to the de Sitter case, where H is a. 
constant.) The time Tc is still a singular point of l|22|l . 
However, it is now a regular singular point since the 
Taylor expansions of the term 



T - Tc{k) 



1 + IbH^fi + 9H^Pt 



and the term 



{r - TcCk)) 



2t(1 + t) 



^Note that to distinguish between functions in rj as opposed 
to functions in r, we will write a tilde on top of functions in r, 
so a = a{rj), but a = a(T) etc. and especially ipj^ = 'Pj^iv) but 

~ (''")■ This should not be confused with the tilde on k 
which was introduced in |^ to denote coordinates in which the 
Fourier modes decouple. 
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indeed exist. For later convenience, we shift the r- 
variable once more so that a; := r + 1; creation time 
for aU modes is now Xc = 0: 



2x{x-l) 



1 + 6H^(3 + 9H^(3x 
AH^(3{x - 1) 



(23) 

Having transformed the creation time into a regular 
singular point, we are now ready to solve the mode 
equation through the Frobenius method. 

4.2 The Frobenius solutions 

Full details on the Frobenius method are given in 
Appendix 1X1 Using the Frobenius ansatz (see H53|) in 
the Appendix) we obtain from H23|l : 



a„(n + r)(3 — 71 — r) • x" ^ 

n=0 



71=0 

oo 

E' 

n=0 



an{n + r){n + r - 1/2) • 
1 \ „_|_i 



+ 3/2 -a; 



n=0 



9/4- > ■«„ -x^^^ =0 



(24) 



The first line of 1)24(1 yields the indicial equation (|52|l . 
which here has the solutions ri — 2> and r2 = 0. 
The first coefficient oo is an arbitrary normalization 
constant and we set ao = 1. Comparing the terms of 
equal powers in x, we then find 



ai = 



02 



r(2r - 1) 



2(r + l)(r-2)' 
1 



4(2 + r)(l - r) 

1 „ r(2r - l)(2r' + 1) 
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From on, all terms in (|24() contribute to the coeffi- 
cient, thus we can find a recursion formula for a„ for 
n > 3 as a function of a„_i, a„_2 and a„_3 as well as 
r: 



a„(r) 



1 



[n + r)(n + r — 3) 



(25) 
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1 



+ -a„_i(n + r - l)(2n + 2r - 3) 



The first solution to l(23() is simply given by inserting 
r = ri = 3 in the above expressions. 



Since ri — r2 = iV = 3 is a positive integer, the so- 
lution for r2 = is of the form given by (|54|l in the 
Appendix. The constant A can be calculated from 03 
as ^ 

^= l™o'''''3W = ^]p;g-«o. (26) 

The coefficients c„(r = 0) for n = 0, 1,2, . 
lated to the a„ in H25|) by 



are re- 



c„(r2 = 0) 



4- ii^ - ^2)a„(r)) 
ar 



(27) 



r— r2— 



The second solution to is then obtained by set- 
ting r = r2 = 0. (Equivalently, A and the coeffi- 
cients Co can be calculated directly by inserting (|5^ 
into H24|) . see Appendix^) The two independent 
solutions can be combined to form the exact general 
solution of 



4>~^{x) = Ci4>J^Jx) + C20fe,2 



(28) 



It is clear from Frobenius theory that the power se- 
ries solutions converge for all a; < 1,, i.e. all r < 0. 
Here, Ci and C2 are complex constants whose choice 
corresponds to the choice of vacuum. Note that our 
solutions H28|l could be rewritten in terms of the con- 
formal time 77 and would then of course solve equation 
®. 

4.3 Wronskian condition 



After the transformation rj 
dition 1)1 3|l takes the form: 



T the Wronskian con- 



i>(T, k) 



V'(r, k) 



(29) 



Note that while the solution space to the mode equa- 
tion has two complex and therefore four real dimen- 
sions, the Wronskian condition reduces the dimen- 
sionality to three real dimensions (since only the imag- 
inary part of H29|) is nontrivial). 

For example, in the case of de Sitter space and 
when splitting the constants Ci and C2 in H28() into 
their real and imaginary parts, Cj = Rj + iKj, j = 
1,2, the Wronskian condition at the creation time 
Tc — —1 can be expressed as: 



R1K2 - K1R 2 = 



12H 



(30) 



4.4 The generic case 

We have here used a suitable transformation of the 
time variable in order to turn the creation time into 
a regular singular point so that then the Frobenius 
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method yielded the exact solutions to the mode equa- 
tion in de Sitter space. In fact, this method, us- 
ing in all cases the same transformation from con- 
formal time to the r variable, can also be applied 
to more general FRW spacetimes. In particular, we 
have applied the method to the case of general power- 
law backgrounds. The exact solutions are then only 
slightly more difficult to calculate, see Appendix [O 
Namely, while in the de Sitter case we obtained a 
recursion formula for the coefficients func- 
tion of a limited number of predecessors (namely three), 
in a general power-law background each coefficient 
generally depends on all others that precede it. In- 
terestingly, the roots of the indicial equation in the 
general power-law case turn out to be the same as for 
the de Sitter case, namely ri = 3,r2 = 0. Therefore 
the behaviour of a mode close to its creation time is 
qualitatively the same in both the de Sitter and the 
general power-law case. 

Further, it is often useful to employ field redefini- 
tions, for example in order to eliminate friction terms 
in the mode equation. Our method remains applica- 
ble under such field redefinitions. In particular, as 
shown in 0, the field redefinition V'fc(^) — V^4'kiv) 
eliminates the friction term in 0, so that we arrive 
at an equation of the form + uj'^{T])ipj^ — 0. Sim- 
ilarly, it is possible to eliminate the friction term in 
our mode equation l|23() in which the creation time is 
a regular singular point. Thereby, the creation time 
remains a regular singular point. The corresponding 
field redefinition and exact solutions are presented in 
Appendix IbI 

Note that formulating the mode equation in fric- 
tionless form allows us not only to obtain the ex- 
act solutions through the Frobenius method but also 
to obtain a second set of solutions using the WKB 
method in the adiabatic regime. This will provide us 
with two bases of the solution space. These are of 
course related by a Bogolyubov transformations and 
we will calculate it Sec|Sl 



tions towards the creation time reveals that the situ- 
ation is more complex. 

Indeed, the values of 4'j.{ii) and its derivative with 
respect to conformal time, at creation time are given 
by: 



Co 



(31) 
(32) 



This shows that prescribing (pj^iji,.) and 'P'^iVc) will 
fix only C2 but not Ci. In fact, at the creation 
time, it is also not possible to specify a solution by 
specifying higher order derivatives of since, as one 
readily verifies, the behaviour of the exact solutions 
shows that all these higher derivatives are necessar- 
ily divergent. For example, </>'^!(?? = f?c) diverges as 
Roughly speaking, as we ap- 



ex 



(plog(-/3feVo2)-l-l) ■ 

proach creation time, the general solution 1)28(1 loses 
its dependence on Ci. 

While the missing Ci has two real dimensions, one 
of these is fixed by the Wronskian condition. Within 
the remaining one-parameter family of solutions, a 
solution can be picked by choosing a phase condition, 
for example, by imposing that the solution be real at 
a given time r^aux ^ ?yc, i-e., Im {(piiriaux)) = 0. 

Thus, any mode function can be specified by giv- 
ing the mode function's value at creation time rjc (two 
real parameters) and demanding that the solution be 
real at an auxiliary time rjaux (one parameter). The 
Wronskian condition constrains the fourth parame- 
ter. 



6 Comparison with the results 
of Easther et al. 

In earlier work, see [7], the irregular singular point of 
the mode equation © was dealt with by truncating a 
series expansion of the coefficients of the mode equa- 
tion. For the resulting approximate mode equation^ 



5 Initial conditions 

Having derived the space of exact solutions to the 
exact mode equation, it becomes crucial to identify 
the physical solution in that space, since this then 
implicitly identifies the vacuum state. 

For this purpose, given that the mode equation is 
of second order, it would appear to suffice to spec- 
ify suitable initial conditions at the mode's creation 
time, such as the physical solution's value and deriva- 
tive. Interestingly, the behaviour of the exact solu- 



(33) 

?7 = ?7c is a regular singular point. Here, jz/ = 5^^^ (H^ 
6(3 H^). The solution space of the approximate mode 



®In |7| a different notation for this equation was chosen; 
apart from replacing the field t/ik by Uf^ = t/ife/a^, also the 
variable y was introduced by setting rj = r]c{l — y), and all 
expressions were re-written in terms of y. 
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equation H33|) was found to be spanned by 



j h \/r] - rjcj cxp{-2^^/£/{'l] - rjc)) 

(34) 

and its complex conjugate, so that the general solu- 
tion is given by the linear combination: 



k^approx 



(35) 



Since we have found in Sec 14. 21 the exact solutions 
(I28|l to the exact mode equation (jS]), we can now 
compare them with the solutions to the approximate 
mode equation that were obtained in f^. 

To compare the behaviour of the exact and the 
approximate solution close to the creation time, we 
take the first and second rj derivatives of F{ri) and 
evaluate at creation time, which gives 



C^^i+^2), (36) 



^k,approx \ '' ' 2 
Aappro.(V ^ ^c) - ^^/^(^l+^2). (37) 

We encounter the same behaviour we observed in the 
case of the exact solution (rj) : the function and its 
first derivative become proportional to another for 
rj ^ r]c and the proportionality factor is the same 



as m (ETJ and namely 2^/ 



2^H * 

Thus, also the solutions to the approximate mode 
equation show that posing initial conditions on (p and 
its derivative at rjc does not suffice to specify a solu- 
tion. 

Further, the second derivative of the solutions to 
the approximate mode equation diverges for 77 — > rye, 
namely as F"{rf) cx In comparison, the sec- 

ond ry-derivative of the exact solution 4>i^{ri) at cre- 
ation time diverges as oc (j)\og{— pk"^ H'^r]^) + 1^ — 

^^2" ^ri-ri + ^(1) (using the series expansion of the 
plog function). Thus, the solutions to the approxi- 
mate mode equation also properly reproduce the lead- 
ing divergence for rj ^ t]c. 

We note that the next to leading orders are dif- 
ferent: the approximate equation has a regular sin- 
gular point at creation time which meant that the 
Frobenius method could be used to derive two power 
series solutions. The roots of the indicial equation 
are (0,3/2), unlike the roots (0,3) of the exact mode 
equation. This means that the solutions to the ap- 
proximate mode equation are of the form: 0^, ^^pp^^^ ^ - 

Eaniv- VcT and (bj^^.pp,,,^^ = Ec«(?7 - Vc)"+^^^ 
The absence of the logarithm in these functions im- 
plies that their higher derivatives exhibit a different 



divergent behaviour than that of the solutions to the 
exact mode equation. 

Related to this, it is clear that the approximate 
mode equation does not conserve the Wronskian of 
the exact mode equation. Instead, it conserves the 
expression — |'^2p- Indeed, the Wronskian of 

the exact mode equation, H13() . for the solutions to 
the approximate mode equation is given by 

= 2Kr7,fc)=0^'/'\/^^(l'^i|' - l'^2p) , (38) 

where 1/(77, given by (UJ, has a nontrivial time de- 
pendence. Using again the series expansion around 
rjc, we find that: 



1 



1 



pl0g(-/3fc2i/2^2) I 



Thus, at r]c, it is possible to enforce the exact Wron- 
skian condition on the solutions to the approximate 
mode equation, though the Wronskian will be con- 
served only to first order in ^/r] — rjc. This shows that 
if a solution to the approximate mode equation is to 
be matched to the numerical evolution of the exact 
mode equation, as was done in 0, then this match 
needs to be performed at a time rjaux which is as close 
as possible to the creation time. We note, however, 
that as our exact solutions will show below, there is 
limit to how close r]aux can be chosen to the creation 
time, because of unavoidable numerical instabilities. 

In j7j , a particular choice of initial conditions was 
suggested, based on a formal similarity between cer- 
tain solutions to the approximate equation and the 
solutions that correspond to the Bunch-Davies vac- 
uum in inflation without a cutoff. The argument is 
that, if the "friction" term in the approximate mode 
equation could be ignored, then it would be of 

the form < 



. 4>l+ujl{vi)(l)k = with ujk^ 



liv) 



In the 



< 1 is 



region where the adiabaticity condition 

satisfied, this equation has two approximate solutions 
of the WKB form: 



k,approx 



4^ 



1/4 



X exp 



(±21^^(7^ - 71c)) .(39) 



Even though one cannot actually ignore the "friction" 
term in the approximate mode equation because the 
period close to rjc is not adiabatic, it is suggestive 
that these would-be WKB type functions possess the 
same oscillatory behaviour as Easther et aZ.'s solution 
to their approximate mode equation. 

In formal analogy to choosing the Bunch-Davies 
vacuum solution, Easther et al. therefore choose "^2 = 
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so that their preferred physical solution to their 
approximate mode equation reads^ 



-^1 



+ is^ \/ri — rjc 



'k^approx " ^ \ 2 

X exp(-2i^£/(?7 ~ rjc)). 

The real and imaginary part of "^i are related via 
the Wronskian condition. The remaining arbitrari- 
ness is merely the freedom to choose an overall phase, 
and can be fixed by demanding, for example, that 
^k,approx ^6*^1 auxihary time -qaux- 

Let us now investigate which exact solution to 
the exact mode equation Easther et a/.'s choice corre- 
sponds to. Normally, in the case of second order wave 
equations, in order to match two solutions at any ar- 
bitrary point it is sufficient to set these solutions and 
their first derivatives equal at that point. We would 
obviously like to match Easther et al.^s function to 
one of our exact solutions to the exact equation at 
the creation time. However, such a match-up at cre- 
ation time itself does here not suffice to use Easther 
et al.^s choice of mode function to uniquely pick out 
an exact solution to the exact mode equation. This 
is because, as we saw earlier, at the creation time 
rj = rjc, the value and the first derivative of the so- 
lutions are not independent (but are instead a fixed 
multiple of another). 

We therefore studied matching-up the mode func- 
tion chosen by Easther et al. with an exact solution 
by setting the functions and their derivatives equal 
at times rjm > Vc- We were indeed able to reproduce 
the observation made in "7? that at times consider- 
ably later than rjc the mode with ^2 = shows close 
to adiabatic evolution with small characteristic oscil- 
lations. This match-up is very delicate, however, and 
we can now see the underlying reasons. 

First, the match cannot be improved by choosing 
earlier and earlier match-up times 77^ because the 
mode function's amplitude and derivative loose their 
independence for 77,„ — > 77c, leading to numerical in- 
stability of the match-up. 

Second, the match also cannot be improved by 
choosing late matching times rjm because the simpli- 
fied mode equation obeyed by 0^, ^^pp^g^ (??) is close to 
the precise mode equation only for early rj. Indeed, 
recall that the function chosen by Easther et al. con- 
serves the Wronskian of the exact mode equation only 
very close to 7]^. 

Given that there is no ideal match-up procedure, 
one may alternatively follow a strategy already out- 

^Note that even though this solution resembles that of an 
adiabatic vacuum, it is not because the adiabaticity condition 

with wj. = (^f^ ) is not satisfied. 



lined in Sec[31 the amplitudes of the two solutions 
are matched at creation time, which fixes C2. The 
information about Ci is then obtained by enforcing 
the Wronskian condition and by choosing the overall 
phase such that the mode function is, for example, 
real at some auxiliary time. We observed again a 
notable dependence on the arbitrary auxiliary time 
Vaux' it must be chosen not too late after creation 
time because the approximate mode equation that 
the function of Easther et al. obeys is close to the 
exact mode equation only in the vicinity of rjc, but 
neither must ijaux be chosen too early because close 
to rjc the exact solution is entirely dominated by one 
of its dimensions in the solution space and informa- 
tion about the contribution of the other dimension is 
numerically increasingly difficult to extract. 

7 Comparison with the WKB 
solution 

The starting point of our present analysis was the 
mode equation ijSJ), which is analogous to that of a 
damped harmonic oscillator. We expressed © in 
the r- variable to obtain (|22|) . which we were then 
able to solve by means of the Frobenius method. As 
we show in Appendix^ a suitable field redefinition 
brings H22|l into a frictionless form. 



4' + (i2(^)^- ^0, 



(40) 



see (|37|l . where 



ci^(r) = 



+132 (3t^H^ 



+5f3H^ + 4:T^ + 8t^ + 4:T 



(41) 



Also this mode equation can be solved with the Frobe- 
nius method. We write the general solution in the 
form 

7^^(T)=i?i4,(r)+i?2V;^2(r), (42) 

see (|59ll in Appendix IbI 

The advantage of the frictionless formulation is 

< 1 



that in the adiabatic range of r, i.e, where 

is fulfilled, (|57|) will have a set of solutions of the 
WKB form, namely 



1 



exp ( +i / u}{T')dT' 



(43) 



A general solution to H57|l in the adiabatic regime is 
then given by ^'^(t) = A_I'^(t) -f A+*t(r). Note 
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that the Wronskian condition, when evaluated at Tc = 
— 1, imposes: 



In the adiabatic regime, the exact solutions and the 
WKB solutions must be related through a Bogolyubov 
transformation. The WKB solutions, H43() . suggest a 
special choice of the vacuum state, namely the adi- 
abatic vacuum characterized by A_ = 0. The Bo- 
golyubov transformation will then tell us which com- 
bination of the exact solutions ■f/'^ i (t) and -0^ 2 (''') 
this choice corresponds to, i.e., it will determine the 
coefficients Di,D2 in (|^ . 

For the purpose of plotting the properties of the 
exact solutions, their Frobenius power series expan- 
sion, lfCT|) . needs to be truncated at some finite order 
N. We established that choosing any iV > 30 suffices 
to ensure that the plot of the truncated function is 
valid well into the adiabatic regime. 

In order to find the Bogolyubov transformation 
between the two sets of solutions, it is sufficient to 
match the value and derivative of the adiabatic solu- 
tion to an exact solution at some arbitrary time Tad 
in the adiabatic regime: 

Dli>k^l{rad) + D2i^k,2i'^ad) 

^X^^^{rad)+X+^{rad) 
Dli^l^iTad) + D2^l^{Tad) 

= A_*'^-(r„<i) + A+*^ + (T,rf) 

We verified that the result of the match-up is in- 
deed independent of the match-up time within the 
adiabatic regime. 

Adiabatic vacuum 

We are particularly interested in that mode function 
V'fe (''') which corresponds to the adiabatic vacuum (to 
zeroth order) in the adiabatic regime. This choice of 
vacuum is defined by A_ = since in the adiabatic 
limit only the modes of positive frequency are present. 
From the match-up of our two sets of solutions we 
then have 

D2 = V/,] • (a+*'+ - i^i^g . 

(Here, the fields are understood to be evaluated at the 
arbitrary time Tad-) The resulting values of Di and 
D2 are several orders of magnitude apart and not ob- 
viously special, in the sense that they do not suggest 



a particular mathematical criterion that would single 
out this initial condition and therefore the vacuum 
in a way that would not rely on the adiabatic expan- 
sion. We noticed however, that numerically, Di and 
D2 appear to be at close to 90 degrees when viewed as 
vectors in the complex plane. In the Outlook we will 
comment on the potential significance for the choice 
of vacuum of the existence of an orthonormal basis 
that is canonical in the solution space. 

Deviations from the adiabatic vacuum 

We can now turn this line of argument around and 
use the knowledge of the exact solutions to bridge the 
gap between the initial behaviour right at the creation 
time and the behaviour of the modes in the adiabatic 
regime. This means that we can precisely link the 
initial behaviour to the one at horizon crossing which 
then determines the size and shape of potentially ob- 
servable effects in the CMB. 

Concretely, starting with a certain linear combi- 
nation of the Frobenius solutions (that may be set by 
quantum gravity) , how close to the adiabatic vacuum 
will this particular solution be during the adiabatic 
phase? To this end, we will derive an expression for 
A_, which measures the deviation from the positive 
frequency adiabatic solution as a function of Di, D2. 
From the match-up (|45|l we find 

+ D^{^l2^t-^t'^L2))- (45) 

So far, A_ depends on the real and imaginary parts 
of Di and D2, i.e. on four real parameters: 

Di = 5i + iLi 
D2 = S2+ iL2 

Only two of these parameters are independent, how- 
ever, due to the Wronskian condition 

^^^-^^^^-(^)"' (46) 

see Appendix IbI and due to the arbitrariness of the 
overall phase, which allows us to choose for example 
Di real. |A_p as a function of these remaining two 
parameters measures the non-adiabaticity of the vac- 
uum in terms of the initial behaviour of the physical 
mode function. Its minimum is zero and denotes the 
adiabatic vacuum. 
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8 Behaviour of physical quanti- 
ties near creation time 

So far, we have parameterized the initial behaviour of 
modes and used our exact solutions to link it to the 
mode's behaviour at horizon crossing. In order to 
better understand what might determine the initial 
behaviour of modes let us now use our exact solu- 
tions to investigate the behaviour of physical quan- 
tities such as the mode's field uncertainties and its 
Hamiltonian close to the mode's creation time. 

8.1 Breakdown of the particle picture 

The Hamiltonian expressed through creation and an- 
nihilation operators was given in H14() . We see that 
the nondiagonal terms of the Hamiltonian, i.e., the 
terms aj,a_j^ and their complex conjugate are pro- 
portional to 



In our case here, we find that at all times rj > rjc 



(^{t], k) = -(/)' 





-kh- (47) 



The condition i^(r/, fc) = can be rewritten as 



= 9 



- ^,{-n,k). 



(48) 



77c, the term 
1 

1 ) . On one 



Let us now recall, see Secl^l that while (/)^ and its 
derivative assume finite values for rj 

virjjk) diverges oc ^plog(— /?^^/a^) 4 

hand, this means that the criterion of 2") for speci- 
fying the vacuum state is here not applicable. On the 
other hand, it is perhaps not surprising that, while 
the modes are created with field fluctuations A^^, of 
finite size, their momentum field fluctuations Att^ are 
divergent at the moment of mode creation itself. 



9 Ground state energy 

Since the Hamiltonian of a mode is quadratic in the 
field TT^ , it is plausible that the divergence of the mo- 
mentum field's fluctuations AtTj, at rjc should imply 
a divergence of the expectation value of the mode's 
Hamiltonian at 77c. To this end, let us consider the 
total ground state energy at some time 77: 



Evaciri) = 



It is clear that at any finite time rj > rjc, this equation 
can be fulfilled, i.e., the Hamiltonian can be made di- 
agonal in the Fock basis, namely by a suitable choice 
of mode function ip^. For the case of the creation 
time 77c itself, however, let us recall from Sec|Slthat all 
mode functions 0^, are proportional to their derivative 
with a fixed proportionality constant, see H31|l and 
(|22Jl. Therefore, (gHIl cannot be fulfilled at the cre- 
ation time 7/c. This means that whatever mechanism 
determines the initial behaviour of modes cannot be 
described in terms of a Fock basis and a correspond- 
ing particle picture. 

8.2 Initial field fiuctuations 

It has been proposed, see , that Planck scale physics 
might imply that when modes are created they ap- 
pear in a state that minimizes the product of the 
fluctuations of the mode's fleld conju- 
gate momentum field tt^. (77) so that the product of the 
uncertainties reads: 




(50) 



Recall that the integration ranges over only those 
modes whose wavelength at time rj exceeds the cutoff 
length. At a given mode's creation time its operators 
ar and at , first enter the total Hamiltonian, which 

k 

means that the mode then first contributes to both 
ground state and dynamical energy. In particular, in 
an expanding universe, the total Hamiltonian contin- 
ually picks up new ground state energy as it picks up 
new modes. 

After inserting fi(r], k) and 1^(7;, k) into (|15|l and 
making use of the fact that we know the exact be- 
haviour of the solutions to the mode equa- 
tion, wc now indeed find a divergence of each mode's 
contribution to the total ground state energy at the 
mode's creation time. Concretely, straightforward 
calculation shows that each mode's Hamiltonian con- 



(49) 
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tains three types of divergent terms: 

cx In (^plog(-/3fcVa2) + 1^ 
cx 1/ (plog(-/3fcVa2) + i) 
cx l/(plog(-/3fcVa2) + i)' 

We obtain a simplified representation of the energetic 
behaviour of each mode by changing variables from 
the pair (j], k) to the new pair (r, fc), and taking into 
account that the integration measure in H5(J|I trans- 
forms as 

T 

In these variables, each mode's ground sate energy 
contains only one divergent term, which behaves as 
1/(t + 1)^ for T ^ —1. We can now consider the 
behaviour of the entity that is of physical signifi- 
cance, namely the integrated ground state energy of 
all modes. 

Here, we notice that the accumulation of ground 
state energy through mode creation is offset to some 
extent by the fact that each individual mode's contri- 
bution to the total ground state energy p^^^ fc(^); i-e-, 
roughly speaking its hioj,/2, continually diminishes as 
its wavelength expands, i.e., as its w^. decreases. 

Whether or not there is a net generating of ground 
state energy through the expansion depends crucially 
on the scale factor function 0(77) of the background 
spacetime. In the case of the background spacetime 
that we have considered so far, namely de Sitter space, 
we know that this calculation can only yield a con- 
stant result, since in its mode equation the depen- 
dence on k can be eliminated, see (|22ll . 

It should be very interesting, however, to carry 
out the integration of the total ground state energy, 
at least numerically, for the case of a more general 
background spacetime, such as the case of power-law 
inflation. As already mentioned, we have calculated 
the exact solutions to the exact mode equation in the 
power-law case, see Appendix O The balance of the 
continual creation of ground state energy and its con- 
tinual dilution through expansion is nontrivial in any 
model of quantum field theory in the sub-Planckian 
regime, including those of generalized dispersion re- 
lations. This is of interest because if the vacuum en- 
ergy generation is not fully offset by its dilution then 
this would imply a vacuum instability and therefore 
a potential starting mechanism for inflation. 



10 Summary and outlook 

We considered quantum field theory in the sub-Planck- 
ian regime, by which we mean the regime of length 
scales larger than but close to the Planck (or other 
cutoff-) length. In the case of an expanding back- 
ground spacetime the independent degrees of free- 
dom of the QFT, i.e., the field oscillators given by 
the comoving modes, are of particular interest. This 
is because as the comoving modes' proper wavelength 
increases, new comoving modes must continually en- 
ter the QFT description of the sub-Planckian regime. 
A key problem is to identify in which state the co- 
moving modes first enter the QFT description of the 
sub-Planckian regime and how they evolve through 
the sub-Planckian regime into the regime where or- 
dinary QFT holds. 

To this end, we here followed up on a concrete 
model for quantum field theory in the sub-Planckian 
regime. In this model a type of natural UV cutoff is 
implemented that has been motivated from general 
quantum gravity arguments, namely the presence of 
a finite minimum uncertainty in positions. Previous 
work had found that in this model the initial phase 
of a mode's evolution is described by a mode equa- 
tion with an intriguing but difficult to handle singular 
point at each mode's starting time. 

Here, we improved on the existing numerical and 
semi-analytical solutions by calculating the set of ex- 
act solutions to the precise mode equation for the 
cases of de Sitter space and power-law inflation. In 
both cases the initial singularity yielded for the roots 
of the indicial equation the values (0, 3). This showed 
that the qualitative behaviour that can be read off 
from of our exact solutions around the initial singu- 
larity is not dependent on the precise dynamics of the 
scale factor. This in turn meant that we were able to 
focus most of our calculations on the simpler case of 
de Sitter expansion. 

Having found the exact solution space, we studied 
possibilities for identifying the physical mode solution 
within that solution space and therefore for identify- 
ing the initial vacuum. The vacuum determines the 
mode's late time behaviour at horizon crossing and 
therefore the type and magnitude of potentially ob- 
servable effects in the CMB. 

Clearly, within any model for QFT in the sub- 
Planckian regime, it is nontrivial to find a reliable 
method for identifying the physical mode solution. In 
particular, given any natural UV cutoff it is of course 
no longer possible to identify a mode's vacuum as 
having started out essentially as the Minkowski vac- 
uum on the basis that the mode would have had ar- 
bitrarily short wavelength in the distant enough past. 
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Within the model of QFT in the sub-Planckian 
regime that we studied here, we found that indeed 
also any approach that is based on instantaneous 
Hamiltonian diagonalization in a Fock basis must run 
into difficulties. The reason is that, as we were able 
to establish using the exact solutions, each mode's 
Hamiltonian enters the total Hamiltonian such that 
it is at the mode's creation time not diagonalizable 
by means of any Fock representation. 

Similarly, the properties of the exact solutions to 
the mode equation showed that the criterion for iden- 
tifying the physical mode function by minimizing the 
field uncertainty relation at mode creation time can- 
not be applied here unchanged: we found that the 
field modes themselves are created with finite fluctu- 
ations. The canonically conjugate momentum field's 
fiuctuations are, however, divergent at creation time 
itself, which is plausible given that the momentum 
field generates changes in the field. 

By studying the behaviour of the exact solutions 
close to the modes' creation time we traced the un- 
derlying mathematical reason for why it is difficult to 
give physical criteria that could reliably single out the 
physical mode solution. Namely, we found that even 
though the mode equation is second order in time, 
at creation time a mode solution cannot be specified 
as usual by giving its amplitude and derivative. This 
is because, at the creation time, the amplitude and 
derivative of all solutions are proportional with the 
same proportionality constant. 

This shows that in any model for QFT in the 
sub-Planckian regime the physical criteria for spec- 
ifying the mode function and vacuum may need to 
be adapted to unexpected mathematical behaviour. 
Here, for example, it is possible to specify every choice 
of mode solution by specifying its amplitude at cre- 
ation time, by using the Wronskian condition and by 
also using the freedom of overall phase to set the am- 
plitude real at some arbitrary time other than the 
mode's creation time. 

Of course it should be most interesting to try to 
calculate directly from candidate quantum gravity 
theories, such as loop quantum gravity and string the- 
ory, the mode solution, i.e., to find the state in which 
these theories predict new comoving modes to enter 
the sub-Planckian regime in which the framework of 
QFT is valid. It would be most satisfying then if 
the fixing of the initial state of the new modes could 
then be rephrased within QFT in terms of a suitable 
boundary action along the lines of |26[ l?7| . Within 
the model of QFT in the sub-Planckian regime that 
we studied here, boundary terms have been consid- 
ered in |2H| and those results are likely to be useful 
for this purpose. 



Independent of which initial condition is to be 
chosen, our exact solutions can be used to straight- 
forwardly calculate from any initial condition the be- 
haviour of the corresponding mode function at hori- 
zon crossing. In other words, the explicit solution 
space provides an explicit bridge between the non- 
trivial initial conditions that may be set by Planck 
scale physics and the predictions for their potential 
impact in the CMB. 

We close with a gedanken experiment that indi- 
cates that Planck scale physics may not necessarily 
set initial conditions for modes that enter the sub- 
Planckian regime during expansion, as this may vi- 
olate unitarity: consider a hypothetical background 
spacetime that repeatedly contracts and expands, so 
that comoving modes repeatedly enter and leave the 
sub-Planckian regime of validity of QFT. In a con- 
tracting phase, when a comoving mode's wavelength 
drops below the cutoff length, the field operators of 
that mode drop out of the QFT's Hamiltonian. This 
means that the further evolution of that mode is 
frozen, at least as far as the QFT is concerned. What- 
ever excitation or particles that mode might have had 
are then unaccessible within the framework of QFT 
because the quantum field simply no longer contains 
operators that act nontrivially on the dimensions of 
that mode's Hilbert space. 

Thus, for all practical purposes any matter in such 
a comoving mode would be disappearing as if behind 
a Planck "horizon" . This would be the case until in 
a subsequent expansion phase the comoving mode's 
wavelength again exceeds the cutoff length. Then also 
the Hamiltonian resumes a nontrivial action on that 
mode's part of the Hilbert space. As far as the frame- 
work of QFT is concerned, if the mode froze while in 
an excited state it will re-emerge in the same excited 
state during expansion: the evolution is unitary and 
no information was lost^. It is possible, of course, 
that the full theory of quantum gravity will show 
that there are operators which act on those frozen 
modes' dimensions in the Hilbert space so that in a 
re-expanding phase, when the modes re-enter the de- 
scription by the framework of QFT, they do so with 
certain fixed initial conditions. Clearly, in this case 
unitarity would be hard to maintain in a cycle of ex- 
pansions and contractions. If, however, Planck scale 
physics were not to enforce a fixed initial condition 
on modes, then the question remains unanswered in 
which state comoving modes are created when they 

*It should be worth investigating if a similar Hilbert space 
mechanism to the one discussed here may freeze degrees of 
freedom that fall into a black hole, to thaw them in the final 
evaporation of the black hole, thus preserving unitarity without 
having the frozen degrees of freedom unduly gravitate. 
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first emerge, during the very first expansion. 

Tliis leads to tlie question whetlier in tfie solution 
space to the mode equation there exists any distin- 
guished or canonical solution that might therefore be 
the preferred physical solution for modes that emerge 
for the first time, i.e., the question is if the mathe- 
matics singles out a preferred vacuum state. This 
amounts to asking if there exists a canonical splitting 
of the solution space that could improve on the use- 
ful but at best approximate splitting of the solution 
space into positive and negative frequency solutions 
during periods of adiabatic evolution. Here we notice 
that, cquivalcntly, the question is whether there is a 
canonical split of the solution space into what in an 
adiabatic phase would be the sine and cosine basis so- 
lutions to the mode equation (from which one would 
then obtain positive and negative frequency solutions 
straightforwardly). Because of the peculiar singular- 
ity at the mode creation time this is indeed possible: 
recall that there exists a distinguished dimension of 
the solution space, consisting of the real- valued mode 
solutions that vanish at creation time. This induces 
a canonical ON basis in the two-dimensional solution 
space that could be viewed as the generalizations of 
the sine and cosine functions. It will be difficult to or- 
thonormalize these functions in practice because the 
exact solution's power series is needed in an integra- 
tion over all time while the power series is slow to con- 
verge at late times. Nevertheless, it will be very inter- 
esting to investigate how close the so-defined positive 
frequency mode solution would be to the adiabatic 
vacuum in the adiabatic phase. 
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A The Frobenius method 

Consider a differential equation of the form 

■ y"{x) +x-[x- p{x)] ■ y'{x) + [x^ ■ q{x)\ ■ y{x) = 0, 

(51) 

where x = is a regular singular point, i.e., the coef- 
ficients have convergent Taylor series expansions 

oo 

X-p{x) = y^PnX" 

n=0 
oo 

and x^ ■ q{x) = QuX^ 

in some interval |a;| < p. The Frobenius method 
states that then there exists a so-called indicial equa- 
tion: 

r • (r - 1) + po • r + go = (52) 

Let ri and r2 be the two roots of this equation, with 
T^i > ''2 , if T^i and r2 are real. Then (|51|l has either on 
the interval — p<a;<0oron0<x<p one solution 
of the form 



oo 

yi(x) = |x|'-^5]a„(ri).a;" 

n=0 



(53) 



In the case where the two roots ri and r2 of H52|) 
differ by an integer (so that ri — r2 ~ N) , the second 
solution takes the form 

oo 

y2{x) = A-yi{x)- In \x\ + \xp ^ c„(r2) • x". (54) 

n=0 

The constant A as well as the a„(ri), c„(r2) can be 
determined by inserting each of the solutions into the 
equation H51l) ^. The solutions converge at least for 
|a;| < p and define a function which is analytic around 
a; = 0. 

^Once having found the coefficients a„ of the first solution, 
A and the c„ can also be calculated from them; see Sec l4.2l and 
Appendix IbI 



B Exact solution to the friction- 
less mode equation 

As shown in 5 , the field redefinition ipj^ (rj) — \/v(f)j^ {i]) 
allows to eliminate the friction term in , so that we 
arrive at an equation of the form + uj'j^{ri)ipj^ — 
with 



/i — 6 
-3 













v' 




< 
















■)-( 


.2v 


\ 



The Wronskian condition (|13|l also simplifies to 

Here we determine the corresponding field redefini- 
tion for the case of the mode equation formulated in 
the T- variable, namely equation (|22|1 . 

Making the ansatz 4>j,{t) = F{T)'ipf^{T) and insert- 
ing into equation H22|l . we find F{t) = ^57^. Recall 
that the range for the variable r is — 1 . . .0. There- 
fore, for later convenience, we take a factor of \/I out 
of F{t) and write the field redefinition in r as 

After the transformation l|55|l , the equation of motion 
for the new field ipj^ir) is tp'^ + ujjXt)%1)j, = with 



W -3 



+3- ' 



- 6 
a / \ a 

5 + T 



- 3- 



ly rf 

a rj' 4t(1 + t) \ D rf 



+ 



5(t + 9) 
16t2(1 + r)' 



(56) 



Substituting a; := t + 1, the mode equation we need 
to solve reads (note that the creation time now cor- 
responds to Xc = 0): 



1 



16/3i7%2(2._i)2 ,36/?ff2a; 
12 ^ H^x^ -27/3 H^x'^ -24/3 H^x 
+32f3H^ ^Ax^ -Ax^ ]iJk 



(57) 







(The prime now denotes d/dr.) Since this equation 
is of the form with p{x) = and the desired be- 
haviour for x^-q{x), we can use the Frobenius method 
to solve it. 
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Application of the Probenius method 

Our ansatz is ipf.ix) = 2;^^^q6„(s) • x", which, 
sorted by powers of a;, leads to 



^ hn [[n + s){n + s - 1) - 2] ; 



sors^° and s: 
bn{s) = 



(n + s)(7i + s - 1) - 2 



n=0 



oo 

+ 

oo 



-2(n + s)(7i + s - 1) + 2 



2(n + 5-l)(n + s-2)-^)6„_i 



{n + s~ 2){n + s - 3) 



27 



1 



27 1 



16 4H'^l3 



K- 



„n+2 



( 1 



Tl = 



9 



On-3 



(58) 



The indicial equation in this case, found from the 
coefficient of , reads 

s(s- l)-2 = 0, 

with the two roots si = 2 and s-i — —1. Note that 
again we find an integer difference between the two 
indicial roots, s\ — si = N = Z. Therefore the first 
solution is of the form (|53|l with s\ =2, while the 
second solution is of the form H54(l with S2 ~ —1. We 
set the arbitrary normalization constant 5o = 1. 
We further find that 



To determine the 6„ for the first solution, simply in- 
sert s = si = 2 in the expressions above. 
The second solution is of the form 

oo 

with Si = —1. Using our results for the first solution 
Vi^i, we find (^ = 3) 



B = lim (s + 1) • bN{s) 



1 



and 



dn{s 



bi 
b2 



4s(s - 1) - 3 
2(s(s + 1) - 2)' 
1 

(s + l)(s + 2)-2 

'(4s(s-l)-3)(4s(s + l)-3) 
4(s(s + l)-2) 

X 27 1 

-s{s-l)- 



(Gs + l)6„(s)) 
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b. = 



1 



(s + 2)(s + 3)-2 



Again, the general solution to (|57|l is given by a linear 
combination, 

i^^^{x)^Dii^^^,ix)+D2i^^,^2ix), (59) 

where Di and Z32 are complex constants. 

Wronskian condition 

To find the corresponding Wronskian condition, we 
plug the field redefinition into (f^ . This leads to 



2(2 + s)(s + l)--j b2 

f , 1^ 27 1 



bi 



(60) 



Splitting the constants Di and D2 into their real and 
imaginary parts. 



For n > 4, we can establish a recursion formula for 
each coefficient 6„ as a function of its four predeces- 



D2 = S2+iL2, 



and 



^"Note that here the regression depth is greater by one than 
in the case of the solution to 1221 . 
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imposing H6U|) at the creation time Xc = leads to 
the condition: 



S1L2 ~ L1S2 — 



12H 



Xc = 0: 



(61) 



G'(g,?7o,0) 



C Exact solution for the gen- 
eral power- law case 

The scale factor in a general power-law background 
reads 



fc 770 e 



00 ^ 

V — 



m— 3 



so that the inverse transformation rj 
(tTC|l ) is given by 



re 



2 \ 2g 



(62) 

T (compare 
(63) 



It can be shown that the exact mode equation 
written for the power-law case still satisfies the nec- 
essary conditions for applying the Frobenius method. 
That is, the creation time Xc = is a regular singular 
point of the explicit transformed mode equation 



Proceeding exactly as before (see Sec0]and Appendix 
, we find the same roots of the indicial equation as 
in the de Sitter case, namely ui = 3, U2 = 0. There- 
fore, the structure of the solutions (i.e., one a mere 
polynomial, (|53(l . the other containing a logarithmic 
term, (I54|) ') is the same in de Sitter case as well as in 
all power-law backgrounds. 

The explicit expressions for the /„((?, 770, u), n = 
1,2,3 found from inserting of the Frobenius ansatz 
are (where we set /q = 1): 



x{x - If 4>l 



(64) 



f2{q,vo,u) = - 



1 

2^ 



x^ - ^x^ + -^a; - 2 



2q 



u{7 - 4m) 
2(l + u)(u-2) 
1 



(u + 2)(w- 1) 

(7-4u)(3/2-2u) 



u \ u — 1 



1 

2^ 



2" / ■^'^ 
+ G{q,T]a;x) (j>f^ 

where we have used the abbreviation 



= 0, 



2{u - 2) 
k'^rjc'^e 



G{q,r]o;x) 



x^ + ^x'^ 



3_ _ 3 

~ 2/ " '4 

3 , kWo 
2 4g2e^-i 



f3{q,Va,u) 



1 



(65) 



I (a;-l)e^-i 
We now work with the Frobenius ansatz 



(3 + u)u 
x\{2 + u) 

{u + 1) 
k'^rjc^e 



-2u 



f2{q,Vo,u) 



i{q,vo;x) = ^fn{q,vo,u)- 



1 



1 

2^ 

/l(«) 



1 

2^ 



k'^rjc e 

4g2 



(66) 



The /n(g,f7o,w) can be determined by inserting H66|l 
into H64|) and sorting by powers of x. Since, how- 
ever, in general the coefficient G{q, rjoi x) need not be 
polynomial, we calculate its Taylor expansion around 



Here we used 77c as an abbreviation; in the power-law 

case rjc — rjo (^/Sk'^e^ . 

In the general power-law background case the co- 
efficient fniq, ijo, u) for n > 4 depends on all the coef- 
ficients preceding it. Explicitly, the recursion formula 
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reads: 

fn{q,vo,u) = - 



1 



{n + u)(n + u — 3) 
x<{n + u—l) 



Using a Taylor expansion of the prefactor, 



X ( - ~2{n + u-2] 



in- 



(n + M — 2) [n + ti — 3 — — 



1/29 



3 ~Pjic^e'^ 

! — 

2 4fl2 



and splitting _Bi,2 into real and imaginary part, Ej = 
Tj + zMj, i — 1,2, the condition we find from evalu- 
ating (|67|l at the creation time Xc = Q reads: 



(- + --3)(i + ^ 



T1M2 - M1T2 



k rjc e 



E 



(m - 1)! 
9™-iG(q,r,o,a:)^ 



m— 4 



fn—r 



The first solution (ui — 3) in the general power-law 
background then reads 

OQ 

4>k,iil^ Vo;x) =x^^ fn{q, r/0, u = 3) x" . 
The second solution has the form 

OO 

where the constant F can be calculated to read 

The coefficients gn{q,r]o,u — 0) of the second power 
series are related to the /„((?, 770, u) by: 



gniq,m,u2 ^ 0) 



-T- (" • fn{q,m,u)) 

du 



- u—u^—O 

The general solution for the power-law background 
case is given by the linear combination of the two 
solutions: 

hi^l^Vo;^) = El cf>f.^^{q,T]o;x) + E2 4>j,2{q,vo; 

El , E2 are complex constants, constrained by the Wron- 
skian condition, which in the power-law case becomes 



2g 



(67) 



i/2g 
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